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Abstract 

The physics of coherent beams of photons carrying axial orbital angular momentum (OAM) 
is well understood and such beams, sometimes known as vortex beams, have found applications 
in optics and microscopy. Recently electron beams carrying very large values of axial OAM have 
been generated. In the absence of coupling to an external electromagnetic field the propagation 
of such vortex electron beams is virtually identical mathematically to that of vortex photon 
beams propagating in a medium with a homogeneous index of refraction. But when coupled to 
an external electromagnetic field the propagation of vortex electron beams is distinctly different 
from photons. Here we use the exact path integral solution to Schrodingers equation to examine 
the time evolution of an electron wave function carrying axial OAM. Interestingly we find that 
the nonzero OAM wave function can be obtained from the zero OAM wave function, in the 
case considered here, simply by multipling it by an appropriate time and position dependent 
prefactor. Hence adding OAM and propagating can in this case be replaced by first propagating 
then adding OAM. Also, the results shown provide an explicit illustration of the fact that the 
gyromagnetic ratio for OAM is unity. We also propose a novel version of the Bohm-Aharonov 
effect using vortex electron beams. 

1 Introduction 

Coherent beams of photons carrying axial orbital angular momentum (OAM), sometimes referred to 
as vortex beams, are well understood, [l] [2] [s] and have various uses in optics and microscopy. [4] 
[5] [6] |7| Recently electron beams carrying very high amounts of axial OAM have been generated [8] 
and the properties of such beams have been studied, [o] [TT] Mathematically the propagation of a 
vortex photon beam in a medium with a homogeneous index of refraction is virtually identical to 
that of a freely propagating vortex electron beam. This is obviously not the case when the electrons 
are propagating in an external electromagnetic field. Here we use the exact path integral solution 
to examine how an electron wave function carrying axial OAM evolves in time. We find that the 



propagation of a wave function carrying nonzero axial 0AM is equivalent to the the propagation 
of a zero 0AM wave function multiplied by an appropriate position and time dependent prefactor. 
Also, the results provide an explicit illustration of the fact the the (non-radiatively corrected) 



gyromagnetic ratio for 0AM is unity as it must be. 11 We will see that from a practical point of 
view this means that the 0AM vector rotates at half the rate of that the electron circulates in a 
magnetic field, i.e., at half the cyclotron or Landau frequency 

The paper is organized as follows Section 2 briefly reviews the derivation of the gyromagnetic 
ratios for orbital and spin angular momentum from the Dirac equation Section 3 discusses the path 
integral solution for the (non-relativistic) propagation of the electron wave function in a magnetic 
field. Section 4 uses the path integral solution to study how a vortex electron beam, actually a 
wave packet, evolves in a magnetic and shows explicitly that the gyromagnetic ratio for 0AM is 
unity. 



2 Dirac to Schrodinger 

For completeness we provide a brief review of the derivation of the Schrodinger equation from the 
Dirac equation which shows explicitly that the (non-radiatively corrected) gyromagnetic ratio for 



orbital angular momentum is unity. 10 
The Dirac equation in SI units is 

{i-ff'D^ - mc) (x, t) = Q (1) 

where ipu is a four-component Dirac spinor and Z?^ — hdf^ — ieA^.Here is the four-vector 
potential and e is the electron charge. The indices fijV, - ■ ■ take the values 0,1,2,3 which correspond 
to the t,x,y,z directions, respectively xq — ct,xi — x,X2 — y,x^ — z. The Einstein summation 
convention wherein repeated indices are summed over their appropriate range is used throughout, 
e.g., Uf,vf = Em=o"m^^^- 

Multiplying Eq (1) by {i^^D^ -I- mc) , and using 

,1 
2 

1 
2 

which follows from {7'', 7"^} — 277''^ where 7^ are the gamma matrices, •q^'^'is the Minkowski metric 



D^D^ - -eha^^'^F^, (2) 



T^^'' = {i/2) [7^,7"^] and F^^, — d^A^ — d^A^^ is the field strength tensor we get 10 

D^'D^ - ^^ehaf^^F^, + m^c^^ i^o {x, t) = (3) 

Consider a constant magnetic field B pointing the in the z direction. Using gauge invariance we 
can write Aq = 0, Ai = —^Bx2 ,A2 — \Bxi^ ^3 = or equivalently Ai = —eij^^xj = —^CijXj. 
Here eijk and e^sare the totally antisymmetric Levi-Civita tensors, e^jfe is +1 (—1) when i,j,k is 
an even(odd) permutation of 1, 2, 3 and is zero otherwise and e^j is +1 (—1) for i, j = 1,2 (2, 1) and 



is zero otherwise 10 Note that diAi — 0. We now have F12 = —F21 = diA2 — 92^i = B. Working 



in the so called "weak field limit", i.e. dropping the A term, gives 
1 



dt^ -d^^] + ^ehB {xid2 - X2di) - eha^^B + rn^c? ) (^,0 = (4) 



2 



In the Dirac basis 



a 




(5) 



10 



where the cr*^ are the Pauli matrices 
and for a slowly moving electron (in the 



In terms of two-component spinors (/) and Xt ''Pd 
Dirac basis) we can set x = and so finally 



(6) 



Here L3 = —ih {xid2 — a;2C^i) is the orbital angular momentum and 53 — ^a^ is the spin angular 



momentum, both in the z direction. More generally 10 

1 



we can write 



(9/-9/ ]-eB- {L + 2S) +m'c^ U(f,t) =0 



(7) 



for a constant B field. Thus we see that the 0AM, L, couples to the magnetic field as B ■ L 
whereas the spin angular momentum, S, couples as 2B ■ S and so the (non-radiatively corrected) 
gyromagnetic ratio for orbital angular momentum — 1 whereas for spin angular momentum 
gs = 2. This difference has the effect that electron helicity, i.e., the spin projected in the direction 
of propagation, remains tangent to the trajectory, i.e, it rotates at the same rate that the electron 
circulates in a magnetic field. We will see below that because = 1 this is not the case for electron 
beams carrying axial 0AM. Note that the values of g^ and gs are a property of the Hamiltonian 
and not of the wave function. The vortex wave function studied below, which carries nonzero axial 
0AM, still couples to the magnetic field with a gL value of unity. 



3 Path Integral Solution for Propagation in a Magnetic Field 

We are interested in 0AM and not spin and so we will drop the spin term in ([t]) and let (x, t) be 
a single component wave function. To reduce to the nonrelativistic case substitute 

0(f,i) = e-™'='*/V(^,i) (8) 

with tp {x,t) slowly varying compared to cxp ^—imc^t/h^ into (j7| and dropping the df.^tp term we 
get the standard Schrodinger equation 

ihdt + —(P + eB ■ L] tljix^t) = (9) 
2m / 

with L — —ihsijkXiXjdk where Xi is the unit vector in the i direction. 

Because ([9]) is linear and first order in the time derivative the solution can be written in the 
form 

^{x,t)= J d^x'K{x,t,x',t')ij{x',t') (10) 

where K {x, t, x' , t') is called the "propagator" and the integral is nominally over all space. The fact 
that ([9]) is first order in time allows the propagator to be written as a path integral 10 12 13 
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I.e., 



K (x, t, x , t') 



5x [t) exp 



dtC (x (t) , dtx (t) , t) 



(11) 



{x',t') 



Here C {x {t) , dtX (t) , t) is the classical Lagrangian corresponding to the quantum Hamiltonian, 
and the integral is over all paths or trajectories which go from x' at time t' to x at time t. The 
Lagrangian corresponding to ^ has the form 



£ {x (t) ,dtx{t),t)^^m {dtx {t)f- eA {x (t) , t) ■ dtX [t) 



(12) 



where A is the vector potential with the magnetic field B = d x A. Using the form for A given 
above we get, for a constant magnetic field in the z direction, 



C {x [t) , dtx (t)) ^ {dtX {t) f + ^eijXidfXj {t) 



(13) 



It should be noted that the Lagrangian in ( 12 ) and ( 13 1 is the full Lagrangian, not the weak field 
approximation . This can be seen simply by calculating the corresponding classical Hamiltonian 

which yields H = (p — eA) /2m. with p = mdtx {t) . 



The solution for the propagator with this Lagrangian is straightforward 12 13 , indeed it's 
given as a problem in Feynman and Hibbs book. 14 Transform to a rotating frame in the xy or 
1, 2 plane by writing 



exp 



eBt 
2to ^ 



cos m 



sm 



f eBt l 



. 2m . 



X2 



(14) 



In terms of the new variables the Lagrangian corresponds to free propagation in the z direction and 
a harmonic oscillator in the Xi, i = 1,2 directions with radian frequency eB /2m. The path integral 
solutions for free propagation and for a harmonic oscillator are well known [12| 13 . Using these 
results and transforming back to the non-rotating coordinates we get 



K{x,t,x',t') = 
with 



3/2 



2mhT) sin [%T] 



exp 



i^z— ^')^ 



cot T] (x, 
mujeijXix'j 



muj 
2 



TO 



(15) 



(16) 



which is the standard cyclotron frequency 13 and T = t — t' . In (15) the combination ujT always 
occurs divided by 2 and so we should expect various aspects of the wave function to evolve at half 
the rate at which the electron circulates in the magnetic field. 
Note that in the limit as cj 



the propagator in ( 15 1 reduces to the free propagator 

3/2 r-„„ „„/\2' 



Kfree {r - ^ ,t ~ t') 



2'Kih {t - t') 



exp 



im (xi 



2h t-V 



(17) 



which is explicitly space and time translation invariant as it should be. 
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Evolution of a Gaussian wave function with and without 
OAM 



The propagator given in ( 15 ) is Gaussian in form and so if we choose a Gaussian for the wave 



function at f = it will remain Gaussian. Also, in this case the integral in ( 10 ) can be evaluated 
analytically. 

First consider propagation perpendicular to the magnetic field. In this case let the initial 
normalized wave function be a Gaussian centered at the origin and propagating in the X2 = y 
direction 

■00 (j, 0) = , , ; exp 



2cr2 2L2 



-py 



(18) 



where we have switched from the Xi notation to the more convenient at this stage x, y, z notation 
with r = cca; + yy + zz. This wave function is roughly a in width in the x and z directions and has 
length L in the y direction. If we specify the values of uj and the radius R of the classical orbit 
of the electron then p = mujR. If we take a and L to be much larger than the nominal de Broglie 
wavelength of 2TTh/p then we expect mininal "diffraction" effects to occur during propagation and 
as shown explicitly below this is exactly the case. This initial wave function has zero OAM about 
it's direction of propagation, the y direction, since 

Ly^o {r, 0) = ih {xd, - zd^) ^Po {r, 0) = (19) 

To generate axial OAM the so called ladder operator approach [l5] is used. Consider an operator 
A with eigenstate \a) such that A |a) = a \a) . We now want to generate a state \a + 1) such that 
A |a + 1) = (a + 1) |a + 1) . To do this we only need to find an operator B such that [A, B] = B since 
then AB \a) = B |a) + BA \a) = (a + 1) B \a) and so the state B \a) = |a + 1) , up to normalization 
and phase factors. Noting that 

[Ly/h, {dx - id:,)] = [i (xdz - zd^) , {d^^ - idz)] {dx - idz) (20) 

it follows that a state with 1 unit of axial OAM, ipi (f, 0) , is given (up to normalization and phase 
factors) by 

V'l (f, 0) - [dx - idz) V-o (r, 0) = ^ {~x + iz) 00 (r, 0) = ^pe^Vo {r, 0) (21) 

Here p — + a,nd 6 increases in the counterclockwise direction when looking in the —y direc- 
tion and is measured from the —x axis. Using the fact that i {xdz — zdx) — —idg we immediately 
see that Lytpi = fi,^i.and so ipi carries one unit of axial OAM. The factor of p, which appears 
automatically, is necessary since at p = (= the y axis in this case) the phase exp [iO] is not defined 
and the wave function must vanish there. 
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Substituting V'o {r, 0) into ( 10 ) and using ( 15 ) gives 
ipo (f, t) = N (fr exp 



where 



2ht 



{z ~ z'f + ^ cot [f ] ({x - x'f + {y~ y'f) 
+ ^ {xy' - yx') 



im 



N exp 

(fir' exp 
-/Vexp 



2;it 4/1 



cot 



y 



Qij^a::' + ayy' + a^z' 



im 9 
z 

2ht 



Ah 



cot 



Lot 

Y 



X (2^)' l3,Py(i, exp 



1 



,/2 



2^/ 2/3, 



N 



/ m \ 
\2iTiht) 



3/2 
' 2 



ttz = 

/3. = 

/3. = 

To propagate tpi we can write 



tmu! 
Y/T 
imuj 

im 

ht 
1 

1 

1 



cot 
cot 



sin[f] / 



2 

y 



2/1 

imu 

^ + yr" 



/7rL2 



imwR 



h 



imio 








cot 




2h 




y 


imuj 






2h 


cot 




y 



^1 (f, t)^N I d\'K (r, r', 0) [d,, - z^,/) Vo (r', 0) 
dVX (f, r*, 0) (-x' + iz') V-o 0) 



T2 



a- 

N 



dx / dVi^ (f, t, r*, 0) exp [A {-x' + zz')] {r^ : 0) 



A=0 



(22) 



(23) 



(24) 



The integral is still Gaussian and can be evaluated as above by letting ax ^ Oix — \ and — >■ a^+iX 
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in (22). Taking the derivative with respect to A and setting A = then yields 



N 

= -exp 


im n 

z + 

2ht 


imu 
cot 

Ah 


Lot 

y 










X ^i2nf 






Pza 


0)exp 


\pxa^ H 


- \Py(^y + \Pza^ 



{-PxOix + ^"00 i) 



(25) 



with UxtPxt ■ -the same as in (23). 

Even though both these analytic solutions can be manipulated into somewhat more convenient 
forms, this is not very illuminating and so we will simply plot these solutions for a set of conditions 
which nicely illlustrate the relevant aspects of their time evolution. On the other hand it is worth- 
while to examine the factor {—PxOix + iPzCtz) to get a better understanding of how it evolves and 
controls the orientation of the 0AM. Substituting from above we find, after some algebra. 



/ (f, t) = -PxUx + iPzOlz = 



[f] 



X + sin [f] y 



sm 



2h 



cos 



1- 



ht 



(26) 



We see that / (r, 0) — —x + at i = 0, as it should, and that it rotates in time in the xy plane 
at a radian frequency of w/2. The origin of this factor obvious. In operator notation, ignoring the 
l/cr^, (21 ) becomes 

l^i) = (-X + *Z)|V^o) (27) 

The time evolution is given by 



e-^"*/'"' (-X + iZ) iV-o) 
(^^-.ut/n + zZ) e+'H*/'') e"'"*/'"' |Vo) 



(28) 



where H = P — eA R j ) /2m is the quantum Hamiltonian corresponding to the Lagrangian 



(13). Note this is the full Hamiltonian, not the weak field approximation. 

The position of the node of ■01 {'r, t) follows from the solution to / (r, i) = 0. At t = this is the 
y axis as shown above. For arbitrary t we have the solution 



' cot 



y 







(29) 



This solution is illustrated in Figure 1 for several values of t. This "nodal line" rotates only by 
TT during one full period, r = 27r/w, of the electron cyclotron orbit and since this factor is the 
origin of the 0AM carried by V'l this shows explicity that the 0AM rotates at half the cyclotron 
frequency, i.e., — I- This also shows that the 0AM is axially oriented only at times t — nr, with 
n = 0, 1, 2, • • • , and its direction switches between being parallel and antiparallel to the direction of 
propagation at each of these times. 
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Note that ipo (^j t) and tpi i^i t) £^re not simply propagating Gaussian envelope functions multi- 
plied by a propagating plane wave factor of the form exp [ip ■ r/h — iEt/H] with |p| constant (but 
rotating at radian frequency uj) and E — \p\ /2m. For both wave functions the de Broglie wave- 
length varies in time. This is to be expected since the coupling to the vector potential contributes 
an extra phase to the wave function of the form —i/h dtA (f) • dtr{t) which varies with position 
in generally an nonlinear fashion . Figures 2 and 3 show slices of the modulus squared and the 
real parts of f/'o and in the xy plane at different positions in the electron orbit. The values 
chosen for (j,L,uj and R are such that the size of the wave packet at i = 0, L in the y direction 
and a in the x direction are both much larger than the wavelength (so that diffraction effects are 
minimal) and R is much larger than L. The actual ratios used for the plots are R — lO'^L, L = 10a 
and a ~ lQ^2Trh/mijj hence the spatial range of the Re [tpo] and Re ["01] plots is about 5 orders of 
magnitude smaller than for the IV'ol^ and 1-0^1 plots so that the phase variation is visible. In Figure 
2 we see that the long axis of the wave function tracks the nodal line and the spatial extent of the 
wave function varies with period r and thus the length and width return, up to diffraction effects 
to their initial values at every t = t, 2t, 3t, • • • . This periodic variation in the spatial extent of 
the wave function can be traced back to the fact that in the rotating frame the Lagrangian is that 
of a harmonic oscillator. The free propagation part of the Langrangian, m (dtx)'^ /2 cause the wave 
function to expand or diffract as it propagates. The harmonic oscillator part, muj'^x^ /2 causes the 
wave function to contract and unless these two effects are precisely balanced the wave function 
will oscillate in size This is exactly analogous to the propagation of a paraxial Gaussian optical 
beam. centered on the z axis and propagating in the z direction in a medium with an index of re- 
fraction of the form n (x, y) = uq — c (a;^ -I- y^) , i.e, a harmonic osciallator potential. In the paraxial 
approximation the propagator for the photon beam has the same Gaussian form as the propagator 
for the harmonic oscillator. The quadratic variation of the index of refraction will case the beam to 
focus or shrink in size as it propagates whereas diffraction effects cause the beam to expand as it 
propagates. If the beam is large, so that the focusing effect dominates, then the beam will shrink 
in size as it propagates. Eventually it reaches a size where the diffraction effect dominates and 
it begins to expand. This process repeats itself causing the beam to oscillate in size with a fixed 



period along its length. 16 These oscillations can be prevented if the size of the beam is fine tuned 



so that the diffraction and focusing effects exactly cancel out. 16 Figure 3 shows the propagation 
of the wave function ^pi carrying a single unit of 0AM. The node in the center of the wave function 
maintains its alignment on the nodal line during each cycle. The spiral form the phase of ipi is 
apparent in the Re [tpi] plots. Clearly the 0AM is rotating at half the cyclotron frequency uj. 



Now consider propagation parallel to the magetic field. In this case we let 

(30) 



■00 {r, 0) , ^ , = exp 



^ + y'^ z^ i 
~2^ " 2L^ + /i^^ 
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Figure 2: Slices in the xy plane of 1-001 and Re [-00] at different positions around the cyclotron orbit 
where -00 is a Gaussian wavepacket carrying axial orbital angular momentum(OAM). The values 
chosen for the width a and length L of the wavepacket, the cyclotron frequency w = eB /m, and 
the radius of the cycloctron orbit R are such that the size of the wave packet at i = (i in the y 
direction and a in the x direction) are much larger than the wavelength so that diffraction effects 
are minimal. All the plots are the same fixed spatial scale with that of the Re [-00] plots being about 
5 orders of magnitude smaller than the |0o| plots so that the phase of the wavepacket is visible. 
At t = 0.5r the wavepacket would be too small to be seen at this fixed spatial scale and so it is 
shown at times t = 0.4r and t = 0.6r instead. 



10 










^0 





Figure 3: Slices in the xy plane of |'0i| and Re [ipi] at different positions around the cyclotron orbit 
where ipi is a Gaussian wavepacket carrying 1 unit axial orbital angular niomentum(OAM) oriented 
in the y direction at i = 0. The values chosen for the width a and length L of the wavepacket, 
the cyclotron frequency lo = eB/m, and the radius of the cycloctron orbit R are the same as in 
Figure 2, i.e., they are such that the size of the wave packet at < = (L in the y direction and 
a in the x direction) are much larger than the wavelength so that diffraction effects are minimal. 
All the plots are the same fixed spatial scale with that of the Re [ipi] plots being about 5 orders of 
magnitude smaller than the \ipi\ plots so that the phase of the wavepacket is visible. Ai t — 0.5r 
the wavepacket would be too small to be seen at this fixed spatial scale and so it is shown at times 
t — 0.4t and t — 0.6r instead. 



11 



and 



ipo {r, t) ^ N d-^r' exp 



- z'f + ^ cot [f] (ix - x'f + iy- y'f 



1 

'20-2 



+ ^ {xy' - yx') 



iVexp 



2ht 



Ah 



cot 







y 





X y d^r'exp 



1 

2^ 



X exp 



™" cot \^] 



4h 



-|3^ 



where TV is the same as in ( 23 1 but now 

/3p = 



1 imu 

1 ini 
1?^ lit 



cot 



im 2 



Ujt 
~2 



(31) 



(32) 



Because depends on a; and y only in the combination = x"^ + y^ it follows that the initial 

Gaussian wave function chosen here does not pick up angular momentum as it propagates along the 
magnetic field. In fact for propagation parallel to the magnetic field the axial 0AM of an eigenstate 
of L2 is conserved. This follows directly from 



[L,,H] =0 



(33) 



where again H = ( P — 6^4 I R ) ) /2r7i and Ai= — ye^Xj. Indeed it can be shown that H 



eB 1 



(X^ + Y^) which obviously yields (33) 



5 Conclusion 

Using the exact path integral solution for the propagator in a constant magnetic field we have 
derived the evolution of a Gaussian wave function and shown explicitly that the (non-radiatively 
corrected) gyromagnetic ratio gi^ for 0AM is unity. This must be the case since is a property 
of the Hamiltonian and not of the wave function. 

The results presented above a novel version of the Aharonov-Bohm effect. [l7j Consider a long 
thin solenoid aligned along the z axis. Outside the solenoid (far from the ends) A varies as 1/ p = 
\l \J x^ + y^ and so B is zero outside. Inside the solenoid A varies as p and so B is constant 
and nonzero. A Gaussian wave function like those considered above carrying nozero 0AM that 
propagates along the z axis has a node on the z axis. In fact wave functions carrying large values of 
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0AM have a very large region around the z axis where the wave function is effectively zero, [s] As 
in the standard Aharonov-Bohm experiment 17 this is a case where there is no overlap between 



the wave function and the magnetic field. The wave function only overlaps with the magnetic 
vector potential. Hence the presence of the solenoid will cause a change in how the wave function 
propagates relative to the no solenoid case. This effect will be predominantly a change in the focus 
position of the wave function. Experimental verification of this would provide yet another example 
of the fact is the fundamental quantity and not E and B. 
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